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Abstract 



Feasibility study of ?7-mesic nuclei production by means 
of WASA-at-COSY and COSY-TOF facilities 



Despite the fact that existence of 77-mesic nuclei in which the 77 meson might be 
bound with the light nucleus by means of the strong interaction was postulated already 
in 1986, it is still not experimentally confirmed. Discovering of this new kind of an 
exotic nuclear matter is very important as it might allow for better understanding of rj 
meson structure and its interaction with nucleons. 

The search of the 77-helium bound states is carried out at the COSY accelerator in the 
Research Center Jiilich in Germany, by means of the WASA detection system. The 
search are conducted with high statistic and high acceptance for the free production 
of the 4 He-77 bound states. It is also considered to search for 77-tritium in quasi free 
reaction which might be realised with COSY-TOF facility. 

In this thesis the results of the Monte Carlo simulations of the 77-helium bound 
states and 77-tritium bound state are presented and discussed. The acceptances of the 
WASA-at-COSY and COSY-TOF detectors for the free and quasi-free ?7-mesic nuclei 
production reactions were determined, respectively. Furthermore acceptances were com- 
pared for three different models of nucleon momentum distribution inside atomic nuclei 
and three different values of width of a considered bound states. In case of COSY-TOF 
detector it was established that the most effcient measurement of quasi-free dd — > p sp (T- 
v)bs — > Ps P dpir~ reaction can be done at beam momentum of pf, eam =3.1GeV/c. 
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1 Introduction 



The new kind of nuclear matter consisting of nucleus bound with r\ mesons via strong 
interaction was postulated by Haider and Liu over twenty years ago pQ. However, till 
now none of experiments confirmed empirically its existence. This exotic form of matter 
called 77-mesic nucleus is schematically presented in Fig. 




Figure 1.1: The scheme of r/-mesic bound state. 

Up to this time ?7-mesic nuclei have been searched via production of r] meson in the 
vicinity of the heavy nuclei. It was considered that due to a huge number of nucleons, 
attraction of r\ meson could be high enough that allows to form a bound state. Nev- 
ertheless those experiments have not brought expected effect. Recent investigations 
indicate that interaction between 77 meson and nucleus is considerably stronger than it 
was predicted earlier. Therefore, according to the current theoretical considerations, it 
is possible that the bound states might be also formed for a light nuclei like helium, 
tritium [21 E] and even deuteron [3j. 

The existence of //-mesic nuclei allows to investigate interaction of the 7] meson and 
the nucleons inside a nuclear matter. Moreover it would provide information about 
N*(1535) resonance [5] and about 77 meson properties in nuclear matter [6j, as well as 
about contribution of the flavour singlet component of the quark-gluon wave function 
of 77 meson [8] . 

The measurement of the 4 He-r/ bound states is carried out with a unique accurance 
by means of WASA detector [D] installed at cooler synchrotron COSY in the Research 
Center Jiilich. The //-mesic nuclei is searched there via studying of excitation function 
for the chosen decay channels of the 4 He-r/ system formed in deuteron-deuteron colli- 
sion [10] • The measurement is performed for the beam momentum varying continously 
around the threshold. The beam ramping technique allows to reduce the systematical 
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Introduction 



uncertainities. The existence of the bound system should manifest itself as a resonance- 
like structure in the excitation curve of eg. dd — > (^He-n)^ — > 3 Hep7r~ reaction below 
the dd — > 4 He-?7 reaction threshold which allows to determine the binding energy and 
the width of such state. 

Formation of the 77-mesic nucleus might be also realized by means of the quasi-free 
reactions. In this case the scan of the energy can be achieved from the Fermi motion 
of nucleons inside the deuteron beam. Measurements of such reaction is available for 
the external COSY-TOF detector [JTj [TJ1 US] where the search of 77-mesic Tritium can 
be carried out by the measurement of the excitation function of the nd — > (T-r])^ — > 
dpit~ reaction around the threshold of the nd — > T-n production [TO] . 

The main aim of this thesis is a determination of geometrical acceptance of WASA 
detector for four reactions in which 77-mesic bound states might be formed via free 
production: 

dd ->■ ( 4 He-r/) fes -)• 3 Hep7r" 
dd — > ( 4 He-r/)fe s — > dppn~ 
pd — > ( 3 He-?])b s — > dpn° — > dp'yy 
pd — > ( 3 He-7])b s — > pppn~ 

and geometrical acceptance of COSY-TOF detector setup for one quasi-free reaction 
of 77-mesic nuclei production: 

dd -)• Ps P {T-i])bs ->■ Pspdpir^ 

In each case Monte Carlo simulations of r/-mesic nucleus production and decay process 
were carried out based on reaction kinematics. The simulations were realized with as- 
sumption that the bound state has a resonance structure given by the Breit-Wigner 
distribution with fixed binding energy B s and a width T. Moreover, for reconstruction 
of events the spectator model was applied. The efficiency of the registration of each 
reactions was analysed with regard to different models describing Fermi momentum dis- 
tributions of nucleons inside deuteron, helium and tritium nuclei. It was also compared 
for different values of B s and T from range predicted by theory [14] . 

This thesis is divided into seven chapters. The second describes indirect and direct 
experimental indications for the existence of the 77-mesic helium. 

The Chapter 3 treats of nucleon momentum distributions inside the light nuclei 
such as 3 He, 4 He and T. The distributions are presented and compared for different 
models. 

Description of the spectator model assumptions and its experimental confirmations 
are presented in Chapter 4. 

The Chapter 5 is devoted to the kinematics of the free and quasi-free reactions in 
which 77-mesic bound states are produced and decay. 
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The WASA-at-COSY and COSY-TOF detection systems emphasising their prop- 
erties useful for the measurement of the 77-mesic nuclei are presented in Chapter 6. 

Simulation results of the bound states in free and quasi-free reactions are described 
in Chapter 7 and in Chapter 8 a summary and conclusions are presented. 



2 Indications for the existence of ?7-mesic helium 



2.1 Indirect 

According to the theoretical considerations, the formation of the ?7-mesic nucleus can 
only take place if the real part of the //-nucleus scattering length is negative (attractive 
nature of the interaction), and the magnitude of the real part is greater than the 
magnitude of the imaginary part [15J: 

\R&{fln— nucleus) | ^* \lTflicin_ nucleus) \ ■ (^-l) 

A wide range of possible values of the r]N scattering lenght a n ^ calculated for hadronic- 
and photoproduction of the rj meson has not exluded the formation of //-nucleus bound 
states for a light nuclei as 3,4 He, T [2], [3] and even for deuteron [I]. Those bound 
states have been searched in many experiments. However, none of them gave empirical 
confirmation of their existence. There are only a signal which might be interpreted as 
an indications of the //-mesic nuclei. 

Experimental observations which might suggest the possibility of the existence of 
the 3 He-r/ bound system were found by SPES-4 [B], SPES-2 (T7J, ANKE [Tg] and 
COSY-11 [H] collaborations. In the experiments cross section of the dp — > 3 He-r] 
reaction was measured. In this reaction with the real r\ meson in the final state, the 
bound state could not be produced and therefore obtained results might be treated as 
indirect indications of the ?7-mesic nuclei only. 




AT D (MeV) 



Figure 2.1: Total cross section for pd — > 3 He-?7 reaction measured for eight different proton 
energies above threshold from 0.2 to 11 MeV [17] . Total cross section rises from 0.25 to 0.40 /xb 
rapidly in the range of 2 MeV. 
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The measurements on the SPES-4 spectrometer were realised using the deuteron 
beam accelerated in SATURNE synchrotron colliding with liquid-hydrogen target, while 
in case of SPES-2 the proton beam and liquid-deuterium target were used. The total 
cross section measured by the SPES-2 collaboration for the eight different beam energies 



are presented in Fig. 2.1 According to [2J, the energy dependence of the cross section 
for the dp — > 3 He-?7 reaction is dominated by the strong interaction between rj and 
3 He originating from the strong r/-nucleon interaction leading to the formation of the 
N*(1535) resonance. The data analysis of the close to threshold measurements of the 
total cross section led to the determination of the 7/ 3 He scattering length. The negative 
sign of a real part of the a^ue and its large value equal a v 3 He = (—2.31 + 2.57z) fm [2] 
suggest a possible existence of the (^He-rj)^, although the condition given by (2.1) is 
not fulfilled. 

An indirect signature of the //-mesic nuclei were also searched at the cooler syn- 
chrotron COSY in Jiilich by means of internal COSY-11 and COSY-ANKE detection 
setups with high precision and high statistics. In the experiments the momentum 
ramping technique of the deuteron beam was used that allows to reduce the systematic 
uncertainties. The beam was accelerated slowly and linearly in time, from excess energy 
of Q=-5.05 MeV up to Q=11.33 MeV in case of ANKE experiment [18], while during 
the COSY-11 experiment [19] the momentum was varied in the range corresponding to 
the excess energy from Q=-10 MeV to Q=9 MeV. Both collaborations performed the 
measurement of the excitation function and differential cross section of the dp — > 3 He-?7 
reaction close to the kinematical threshold. The experimental excitation function 



parametrized with the s-wave formula of scattering length [T8l [T9] is presented in Fig. 2.2 
(left panel). The fit to the COSY-11 data gave the value of the 7/ 3 He scattering length 
equal to a v 3jj e = [±(2.9±0.6) + (3.2±0.4)i] fm [19] . Although this value is in agreement 



with formula (2.1 ), uncertainties of its real and imaginary part are too large to confirm 
the possible formation of ^He-r])^. The real part of scattering lenght of the 7/ 3 He system 
derived by fitting the ANKE data for Q<4MeV equals Re{a n i He ) = (11.6±1.4) fm while 
imaginary part is equal Im{a ri 3 He ) = (—4.1 ± 7.0) fm [IB]. Those large values implies 
the existence of a quasi-bound states very close to the reaction threshold, however the 
derived by ANKE and COSY-11 collaborations real parts of the scattering length are 
not consistent within the quoted errors. 

The differential cross section measurement allows to calculate angular asymmetry pa- 
rameter a, defined as: 

d , da , 
d cos u v d\ I 

The momentum dependence of parameter a can be described correctly with assumption 
that the phase is varying between S and P waves [21] . In another case the discrepancy 
between the experimental data and theoretical description are significant. The momen- 



tum dependence of a parameter is presented in Fig. 2.2 (right). 
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Indications for the existence of eta-mesic helium 




Figure 2.2: (left) Total cross section for the dp — > 3 He-r/ reaction measured with the 
COSY-ANKE (open circles) p] and the COSY-11 facilities (closed circles) [H] and (tri- 
angles) [20]. Scattering length fit to the COSY-ANKE and COSY-11 data is represented with 
dashed and solid lines, respectively, (right) Angular asymmetry parameter a for the experi- 
mental data from COSY-ANKE (full dots) [18] and from COSY-11 (open circles) [H]. The 
dashed and solid lines are fitted (assuming the phase variation between S and P waves) to 
the COSY-11 and COSY-ANKE data, respectively [21]. The fit without the phase variation 
is denoted as the dotted line. The figure is adapted from [10J. 



2.2 Direct 

The first direct experimental indications of a light 77-nucleus bound states were 
observed in the reaction of the r\ photoproduction 7 3 He — > n°pX which was investi- 
gated with the TAPS calorimeter at the electron accelerator facility Mainz Microtron 
(MAMI) [22] . Photons produced in a thin radiator foil and tagged using a special spec- 
trometer hit the target filled with liquid 3 He. There the measurements of the excitation 
functions of the 7r°-proton production for two ranges of the relative angle between 
those particles were carried out. It appeared that a difference between excitation curves 
for opening angles of 170° — 180° and 150° — 170° in the center-of-mass frame revealed 
an enhancement just below the threshold of the 7 3 He — > 3 He-r/ reaction which was 
interpreted as a possible signature of a 3 He-r/ bound state where 7] meson captured by 
one of nucleons inside helium forms an intermediate «Sii(1535) resonance which decays 
into 7T°-p pair. A binding energy and width for the anticipated quasibound r/-mesic 
state in 3 He were deduced from the fit of the Breit-Wigner distribution function [22] to 
the experimental points and equal (—4.4 ±4.2) MeV and (25.6 ±6.1) MeV, respectively. 
Those values are consistent with expectations for 77-mesic nuclei. Above cited excitation 
functions and the difference of them with a Breit-Wigner distribution and background 



fitted to the data are shown in Fig. 2.3 
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Figure 2.3: Excitation functions of the 7r°-proton production for relative angles of 170° — 180° 
(red triangles) and 150° — 170° (black circles) in the 7 3 He center-of-mass sytem are shown 
in the left and center panels. In the right panel the difference between both distributions 
with superimposed line denoting the results of the fit of the Breit-Wigner distribution plus 
background are presented. The figure is adapted from [22]. 



However, due to the low statistics of the measurement the results might be interpreted 
not as an indication of the bound state but rather as a virtual state what is in details 
described in Ref. [23]. The interpretation is still under discussion [23]. Moreover at 
the recent meeting it was shown that the result may be an artefact due to the strong 
influence of the resonances on the shape of the excitation function [25] . 




-0.5 -0.25 0.25 0.5 0.2 0.4 0.6 0.8 1 20 40 60 80 100120 140 160 180 



p L [GeV/c] p [ GeV/c ] e cm p-- [ deg. ] 

Figure 2.4: Experimental results of the COSY-11 collaboration for the dp —> pppir~ reaction: 
(left) Transversal vs. longitudinal momentum distributions of protons. The upper limit for 
the longituidal proton momenta is shown as dashed line and equals p£=0.18GeV/c. (middle) 
Pion momentum distribution in the center of mass system, (right) Relative angle between 
pion and proton direction in the cm. The figure is adapted from [27J. 

The analysis carried out by COSY-11 group [2B1 [271 EHl I2H] give an indication 
for the 3 B.e-r] bound state existence. The search for the ?7-mesic helium was carried 
out using a deuteron beam and internal hydrogen target. The beam momentum was 
ramped around the kinematical threshold for the r\ production in the dp — > 3 He-?7 
reaction and the measurement of the dp — > ppprt^ and dp — > 3 He7r° reactions was 
carried out. In the first case the momentum distribution of the 7r~ (middle panel 
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Indications for the existence of eta-mesic helium 



in Fig. |2.4 ) and the relative angle distribution between pion and proton momentum 
vectors (right panel in Fig. 2.4) were determined after application of apropriate cuts 
on the momentum of the spectator protons (left panel in Fig. 2.4) and the rejection 
of the events corresponding to quasi-free 7r~ production. Obtained results are in 
agreement with theoretical expectations for particles originating from decay of the 
N*(1535) resonance which is created as a result of the absorption of the bound t] 
meson in the neutron inside 3 He. Based on the above results the upper limit of total 
cross section for the dp — > ( 3 He-r])b s — > pppn~ reaction was estimated to the value of 
270 nb. Similarly, investigation of the dp — > 3 He-7r° reaction give only the value of the 
total upper limit of cross section of the dp — > ^He-r/)^ — » 3 He7r° reaction equal to 70 nb. 



3 Nucleon momentum distributions inside d, T, 3 He and 4 He 
nuclei 



Due to the Fermi motion, nucleons inside atomic nuclei are not at rest but move with 
momenta which vary in a broad range. This variation influences kinematics of nuclear 
reactions. 

Fermi momentum distributions of proton and neutron bound inside a deuteron de- 
rived from two different potential models, namely PARIS [30] and CD-Bonn [31] are 
shown in Fig. 



3.1 



The normalized nucleon momentum distributions are calculated by 
means of the Fourier transformation of parametrized deuteron wave functions obtained 
from space representation. Respective parametrization coefficients found in the analitic 
representation of the deuteron wave function for both of above named nucleon-nucleon 
interaction models are given in (30j EU [32] . The momentum distributions deduced from 
Paris and CD-Bonn potentials are peaked at about 40 MeV/c and differ no more than 
5% 132]. 




100 200 

Ppermi [ MeV/C ] 



300 



Figure 3.1: Fermi momentum distribution of nucleons inside the deuteron for PARIS (full 
line) and CD-Bonn (dashed line) potentials. The distributions were normalized to unity in 
the momentum range from to 300 MeV/c. 

In case of three-nucleon bound states like 3 He and T, Fermi momentum distributions 



of nucleons are presented in Fig. 3.2 for three different models. Thick solid line depicts 
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Nucleon momentum distributions inside d, T, 3 He and 4 He nuclei 



proton momentum distribution inside 3 He and neutron momentum distribution inside 



T as given by analytic formula (3.1 ) which results from the fit to the experimental data 



on p ( 3 He, 2p) d and p (T, pn) d reactions 

f(p) = p 2 [exp(-263p 2 ) + 0.177exj9(-69.2p 2 )] (3.1) 

This momentum distribution is in a good agreement with the one calculated with real- 
istic potential in the frame of the model of the composite quark bags [53] . 

The distributions of protons and neutrons momentum inside 3 He and T are also 
estimated based on the AV18 and the CDB-2000 nucleon-nucleon interaction models in 
conjunction with Urbana IX (UIX) and Tucson-Melbourne (TM) three nucleon interac- 



tions (TNI), respectively [351 [36]. They are presented in Fig. 3.2 for protons inside 3 He 
(left) and neutrons inside T (right) and ticked as a dashed and dotted lines. Similar 
results can be obtained for neutron inside helium and proton inside tritium. The dif- 
ference between those two distributions are small and results from different interaction 
Hamiltonian forms defined for above-cited models. 




Figure 3.2: Fermi momentum distribution for protons inside 3 He (left) and neutrons inside 
T (right) given by analytic formula (thick line) and estimated for the AV18 NN (dashed 
line) and the CDB-2000 NN (dotted line). The distributions were normalized to unity in the 
momentum range from to 0.4 GeV/c. 



Estimation based upon the formula (3.1) is consistent with the one derived from AV18 
NN and CDB-2000 NN models with an accurancy better than 9% for protons inside 
3 He and 11% for neutrons inside T. 

For nucleons inside 4 He Fermi momentum distributions predicted by three indepen- 



dent models are shown in Fig. 3.3 The distribution represented by a thick line is 
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calculated from helium wave function derived based on Fermi three parameter charge 



distribution of nucleus [37]. The momentum distribution is described by formula (3.2): 



f(p) = — h exp (-£) , (3.2) 
a ■ o \a ■ c J 

where a = 0.03892719, b = 0.05511, c = 0.7352. Fermi momentum is given in units of 
GeV/c. 

The dashed and dotted lines depict distributions obtained, similarly as in case of 
three nucleon systems, from AV18 and the CDB-2000 potential models with the in- 
clusion of three nucleon interaction contributions [36]. Due to the fact that 4 He is 
symmetrical, proton and neutron momentum distributions are in good approximation 
equal. 




Figure 3.3: Fermi momentum distribution of nucleons inside 4 He given by analytic formula 
(thick solid) and estimated for the AV18 NN (dashed) and the CDB-2000 NN (dotted). The 
distributions were normalized to unity in the momentum range from to 0.5 GeV/c. 

The difference between the distributions derived from AV18 and CDB-2000 models 
and given by analytic formula is significant and equals up to about 40%, and in 
addition the maxima of these distributions are shifted by about 45 MeV/c. The 



discrepancy results from the fact that the formula (3.2) was derived from nucleus 
charge distribution smeared out by the charge distribution of protons, whereas the 
AV18 and the CDB-2000 models allow for the finite size of nucleus charge distributions 
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Nucleon momentum distributions inside d, T, 3 He and 4 He nuclei 



and are related to the momentum of the point like protons in the alpha particle [38] . 
Therefore, as more realistic in further considerations the momentum distributions of 
nucleons inside 4 He calculated from AV18 and the CDB-2000 potentials will be taken 
into account. 

The numerical data describing above cited momentum distributions are given in 
Appendix A. 

The presented distributions estimated based on various models will be used to cal- 
culate the systematical uncertainty in the determination of the excitation curves used 
to search for r/-mesic bound states. Respective results are presented in Chapter 7. 



4 Spectator model 



4.1 The main assumption of the spectator model 

The production and decay of //-mesic bound states investigated in this thesis might be 
schematically depicted as: 

3 He sp p it' 



a) «H -(*«,),.-> \ p, j,,r 
c) ->• p sp nd -i> p sp (T-T]) bs -> p sp d sp p 7r~ 



The a) and b) schemes include free and c) one describes quasi-free production of the 
bound states. Subscript bs denotes 'bound state' whereas sp stands for the 'spectator'. 

'Spectators' are particles which do not take part in reactions [39] HQ] but hit the de- 
tectors with the Fermi momentum transformed into laboratory system. Fermi momen- 
tum distributions of nucleons inside the light nuclei are presented in previous chapter. In 
the framework of the spectator model due to the relatively small binding energy of the 
nuclei, spectators are considered as a real particles registered in the experiments and in 
the analysis it is assumed that they are on their mass-shell during the reaction |39j SO] : 

|P sp | 2 = m%. (4.1) 

The F sp and the four-momentum vector of spectator and the spectator mass, 

respectively. 

In the free reactions the beam and target nuclei collide and form ?7-mesic bound state 
which decays into proton, pion and spectator/spectators which energies and momenta 



are measured in experiment. One of those reactions is schematically shown in Fig. 4.1 



In case of quasi-free reaction presented in Fig. 4.2, the deuteron from the beam is 
considered as a system consisting of proton and neutron moving with the Fermi motions. 
For the deuteron beam we have: 

F d = F b n + F Psp (4.2) 
\F Psp +F b n \ 2 = ml (4.3) 
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Spectator model 
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Figure 4.1: Schematic picture of the dd — > ( 4 He-?7){, s — > 3 Hep7r _ reaction. Red and blue circles 
represent protons and neutrons respectively, whereas tt~ meson is depicted as yellow circle. 
The beam momentum is presented by the dashed arrow. 

where rrid denotes the deuteron mass equal to 1875.6 MeV/c 2 [H] while F Psp and are 
the four-momentum vectors of the proton spectator and the beam neutron, respectively. 
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Figure 4.2: Schematic picture of the quasi-free dd — > p S p(T-v)bs Pspdpir~ reaction. Red 
and blue circles represent protons and neutrons respectively, whereas rj meson is depicted as 
green circle. The Fermi momentum of the nucleons inside the deuteron is presented by the 
dotted arrows and the beam momentum by the dashed one. 

According to the spectator model, proton from the beam does not take part in the 
reaction and is registered as a real particle whereas neutron being off the mass-shell 
hits the target deuteron. From the conservations of momentum and energy and the 
assumption that proton is on its mass-shell we obtain in the deuteron beam rest frame: 



Pn = -P. 



(4.4) 



E* n = m d - E 



Psp' 



(4.5) 
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Based on above-mentioned relationship we can deduce the neutron four-momentum 
vector from the spectator momentum which is fundamental in the analysis of reaction 
kinematics. The close description of free and quasi-free reactions processes is described 
in Chapter 5. 



4.2 Experimental proofs of spectator model 

The validity of spectator model assumptions was confirmed by measurements performed 
by collaborations WASA/PROMICE g2], TRIUMF g3], COSY-TOF @3|, COSY- 
11 [45J and HADES [46]. The WASA/PROMICE collaboration [42] has compared free 
and quasi-free production cross sections for the pp — > ppr\ reaction. As a result it was 
presented that within the statistical errors there is no difference between the total cross 



section of the free and quasi-free process. The experimental data are shown in Fig. 4.3 
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Figure 4.3: Total cross section for the pp —> pprj reaction as a function of the excess energy for 
free (open circles) and quasi-free proton scattering (full circles). Figure is adapted from [17] . 
The data are taken from references [Ml HH EH EH E21 E2 [54] . 

Investigation of pion production at the TRIUMF |43J facility in quasi-free pp — > dn + 
reaction extracted from the pd — > dir + n reaction has proven that the spectator momen- 
tum distribution determined from the experimental data agrees with expectations based 
on spectator model. Moreover, it was shown that the magnitude of the differential cross 
sections for the quasi-free and for the free reactions are consistent on the few per cent 
level. 

The spectator assumption was also confirmed by the COSY-TOF group [S|. The 
momentum distribution of the spectator as well as the shape of the angular distribution 
for the quasi-free np — > ppir~ and pn — > pn reactions have been measured. The exper- 
imental data are consistent with calculations based upon the hypothesis of spectator 
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Spectator model 



model with an accurancy better than 4% up to 150 MeV/c of the Fermi momentum 
and with about 25% up to a momentum of 300 MeV/c. 




Figure 4.4: Proton spectator momentum distribution reconstructed in COSY- 11 experiment 
(points) in comparison with simulation taking into account Fermi momentum distribution 
of nucleons inside deuteron, the acceptance and the efficiency of the detector system (solid 
histogram). The figure is adapted from |45j . 

In case of quasi-free pn — > pnr]' reaction studied at COSY-11 facility [45J it is shown 
that the measured proton spectator momentum distribution is in good agreement with 
the theoretical assumptions of spectator model. The comparison of the experimental 



data and simulation result is shown in Fig. 4.4 



Recently, the validity of the spectator model was proven also by HADES collabo- 
ration jl6] during the measurement of the quasi-free np — > e + e~pX reaction realised 
with a deuteron beam and proton target. The angular distribution of proton spectator 
outgoing from deuteron beam and its momentum distribution agrees with assumptions 
of spectator model up to 200-300 MeV/c. Momentum distributions in deuteron CM 
frame are presented in Fig. 4.5| Experimental dat ticked red points, while 



simulations results based on the spectator model as a black points. 

The above-cited results confirmed the spectator model and thus allow to use it in 
the analysis of the reactions kinematics. 



25 




Figure 4.5: Momentum distribution of proton spectator measured in HADES experiment 
(red points) and calculated based upon the hypothesis of spectator model (blue points) for 
M e+e - < 140MeV/c 2 (left) and for M e+e - > 140MeV/c 2 (right). Picture courtesy of 



5 Kinematics of the 7/-mesic bound states production and de- 
cays 



In this thesis four reactions of free and one of quasi-free 77-mesic bound states production 
are considered 



1. dd — > ( 4 He-r]) 6s — > 3 Hep7r~ 

2. dd — > ( 4 He-?7)ft s — > dppm~ 

3. pd —> ( 3 B.e-r]) bs — > dpn° — > dp'yy 

4. pd —> ( 3 He-r]);, s — > pppn~ 

5. nd (T-r/)^ — > dp7t~ 

In case of free reactions (l)-(4), A B.e-r] and 3 He-?7 bound states are produced in 
deuteron-deuteron and proton-deuteron fusion, respectively [SB]. The mechanism of 
the reactions is presented schematically in the example of the ( 4 He-?7)b s production in 



Fig. 5.1 Describing the kinematics of the reaction following notations will be used: 



P^ = (i? (j b ,pfe)-four-momentum vector of the beam deuteron 
Pj = (rrid, 0)-four-momentum vector of the target deuteron 
^ 3 He = (-E3# e ,p3# e )-four-momentum vector of the outgoing 3 He 
P p = (-E , p ,pp)-four-momentum vector of the outgoing proton 
Ptt- = (E n - , p n - )-four-momentum vector of the outgoing pion 



SCHEME OF REACTION PROCESS , 
IN WHICH ri-MESIC NUCLEUS IS FORMED 




c 

I 3 He, PROTON 

AND MESON %~ 
EMISSION 



CREATION OF 
Ti-UESIC NUCLEUS 



RESONANCE DECAY 
INTO PION AND 
PROTON INSIDE 
NUCLEUS 



ABSORPTION OF i, MESON BY NUCLEON EXCITATION INSIDE 

ONE OF NUCLEON INSIDE THE THENUCLEUS- 

HELIUM N* RESONANCE FORMATION 



Figure 5.1: Reaction process of the (^Be-r/)^ production and decay. 
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According to the scheme shown in Fig. 5.1, the deuteron from the beam hits the 



deuteron in the target with a momentum of p b - The collision may lead to the creation 
of 4 He nucleus bound with the r] meson via strong interaction. The mass of a bound 
state is a sum of r\ and 4 He masses reduced by binding energy (B s ): 

m bs = m v + rrn He - B s . (5.1) 
The //-mesic nucleus moves in laboratory frame with velocity: 

where Sdd is the square of invariant mass of the colliding deuterons: 



s dd = \K + K\ = 2m d U + + Wj • (5-3) 

The rj meson might be absorbed by one of the nucleons inside helium and may propagate 
in the nucleus via consecutive excitation of nucleons to the N*(1525) state [57] until 
the resonance decays into the pion-proton pair outgoing from the nucleus [TDJ [561 EH] • 
Before the decay, it is assumed that N* resonance moves with a Fermi momentum pp 
inside 4 He. From the momentum conservation in the 4 He frame and the assumption of 
spectator model, momentum and energy of 3 He may be expressed as: 

Pm e = ~Pf (5-4) 



Es* He = ^jml He + b7l 2 - (5-5) 

The momentum and energy are transformed into the laboratory frame by means of 
Lorentz transformation: 

Vme = Vme + Pcmlcmilcm/ (icm + l)Pcm * PVte + E We) ( 5 ' 6 ) 

— t 

Ei H e = lcm{Ei He + f3 cm ■ p 3 * He ), (5.7) 



where 7 cm = 1/y 1 - |/3 cm | 2 . 

The angle between outgoing 3 He and the beam direction is given by: 

a ( P 3 He-Pb \ , K t 

Uz He = arccos — . — . (o.c 
\\pme\ ■iPblJ 
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The relative angle between the outgoing nucleon-pion pair is equal to 180° in the N* 
reference frame. In the following the variables in the N* reference frame will be denoted 
by '**'. Both particles move with a momentum \p **_ | which is related to the resonance 
mass: 

m N * = (^s dd + ml He - 2 y /s d ~ d ~^Jm 2 He + \pp\ 2 ^j , (5.9) 



and is given by: 



\{m 2 N *, m 2 _, m 2 ) 



where X(x, y, z) = (x — y — z) 2 — Ayz [59|. 

The pion and proton four-momentum vectors in the laboratory frame are calculated 
using the Lorentz transformation, first from N* to the bound state frame: 

p p V = p*:- + n * 1n *(jn*/(jn* + ■ p*1- + e;;-) (5.11) 

e p %-=jn*(e;;- + p n *-p;;-), (5.12) 

and further to the laboratory frame: 

— * — * 

Pp,*- = Pp,*- + Pcmlcm{lc m / (7cm + I) Pan ' Pp,n~ + E p, V -) ( 5 ' 13 ) 

Ep^- = 7cm(E;, n - + (3cm ■ P p V)> ( 514 ) 



where 77V* = 1/y 1 — \(3n* \ 2 is a velocity of the resonance N* in the bound state frame. 
The angles of outgoing proton and pion in LAB system equals: 

a ( Pp,*- -Pb \ /r 1 r\ 

U V7T - = arccos 7— - — . — r— - . (5.15) 

P \\Pp,n-\-\Pb\) 



The quasi-free reaction kinematics was partially characterized in Chapter 4 by the 
way of the spectator model description. Neutron bound inside the deuteron hits the 
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deuteron target forming the (T-r))b s , while proton does not take part in reaction and is 
considered as a real particle: 

KJ 2 = m 2 Psp . (5.16) 

The square of invariant mass of neutron-deuteron system (s n d) depends on the neutron 
Fermi momentum p* (5.1.3) inside the beam deuteron and taking into account that 
target deuteron is at rest in the laboratory s nd equals: 

s nd = (E n + m d ) 2 - |p„| 2 , (5.17) 

where p n and E n are neutron momentum and energy in laboratory frame which may 
be obtained applying the Lorentz transformation according to the following formulas: 

Pn = p: + M d {id/{id + \% ■ p: + ED (5.i8) 

E n = ld (E* n + p d -p:) : (5.19) 



where (3 d denotes velocity of the beam deuteron in the laboratory frame, and p* denotes 
neutron momentum in the deuteron center of mass, and E* n = m d — \jm 2 v + |p* | 2 . 

The process of T-rj bound state decay is analogous like in case of the free reaction 
which kinematics was described before. The deuteron spectator escapes with the Fermi 
momentum and energy, which in the frame of the bound state is equal to: 

p d * = -p; (5.20) 



E* d = yjm* + \pt\*, (5.21) 
After the transformation into laboratory system we have: 

Pd = P d + P cm' 'J cm' (icm'/dcm' + l)Pcm' • P d * + E*) (5.22) 

E d = lcm ,(E d * + p cm ,-p d *), (5.23) 



where Pcm' = m ^+ E denotes velocity of the center of mass for the quasi-free nd — > (T- 
r])b s dpn~ reaction in the laboratory frame. 

The outgoing proton-pion pair originates from the decay of the resonance created 
via absorption of the r\ meson on a nucleon in the tritium nucleus. The four-momenta 
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Experimental setup 



of those particles are described in resonace frame by equations analogous to (5.11) and 
(5.12), and in laboratory frame by formulas (5.13) and (5.14) while proton and pion 
angles with relation to the beam direction are given by (5.15). 



6 Experimental setup 



6.1 WASA-at-COSY facility 

The search for ^-mesic helium in free production reactions with high statistic and ac- 
ceptance is carried out at the WASA facility [9], an internal detection system installed 
at the cooler synchrotron COSY in the Research Center Jiilich. The WASA detector 



vertical cross section is schematically presented in Fig. 6.1 All setup components and 



the method of measurement are described in detail in reference [9J. Thus, in this 
chapter the experimental technique will be only shortly presented. 



Pellet Target 




Central Detector Forward Detector 



Figure 6.1: Scheme of WASA-at-COSY detection system. Gamma quanta, electrons and 
charged pions being products of mesons decays are registered in the Central Detector. Scat- 
tered projectiles and charged recoil particles like 3 He, deuterons and protons are registered in 
the Forward Detector. The abbreviations of the detectors names are explained in the text. 



In the COSY synchrotron protons and deuterons might be accelerated in the mo- 
mentum range between 0.3 GeV/c and 3.7 GeV/c [9J. The ring can be filled with 
up to 10 11 particles leading to luminosities of 10 31 cm _2 s _1 in case of internal cluster 
target [HD] and 10 32 cm~ 2 s _1 in case of pellet target [5j. Beams are cooled by means of 
electron cooling as well as stochastic cooling at injection and high energies, respectively. 

The internal hydrogen (H 2 ) or deuteron (D 2 ) target of the pellet-type is installed in 
the central part of the WASA-at-COSY detector and its position is marked in Fig. 6.1 as 



a dotted line. The central detector is built around the interaction point and designed for 
measurements of 7r° and t] mesons decay products like photons, electrons and charged 
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pions. The charged particles momenta and reaction vertex are determined by means of 
Mini Drift Chamber (MDC) which covers angles from 24° to 159°. Charged particles are 
here bending in the magnetic field provided by sourrounding Superconducting Solenoid 
(SCS). First their trajectories are reconstructed, and then knowing the magnetic field, 
the momentum vector is reconstructed. For identification of charged particles the AE-p 
and AE-E methods are used based on AE signals in Plastic Scintillator Barrel (PSB). 
The photons, electrons and positrons are registered in Scintillator Electromagnetic 
Calorimeter (SEC) via production of electromagnetic cascades. The calorimeter covers 
polar angle in the range from 20° to 169°. 

The detection and identification of forward scattered projectiles and target-recoil 
particles such as protons, deuterons and He nuclei and also of neutrons and charged 
pions are carried out with the Forward Detector which covers the range of polar an- 
gles from 3° to 17°. It consists of fourteen planes of plastic scintillators forming For- 
ward Window Counter (FWC), Forward Trigger Hodoscope (FTH), Forward Range 
Hodoscope (FRH), Forward Range Interleaving Hodoscope (FRI) and Forward Veto 
Hodoscope (FVH), respectively and proportional counter drift tubes called Forward 
Proportional Chamber (FPC). Particles trajectories are reconstructed from the signals 
registered successively in FWC, FPC, FTH and FVH scintillator modules. Particles 
are identified based on measurement of energy loss in the detection layers of FRH, 
FWC and FTH. The registered energy loss allows to determine their total momentum 
which direction is reconstucted from the measurement of particles tracks by means of 
straw detectors constituting FPC. Respective components of the Forward Detector are 



presented in Fig. |6.1 

The 3 He-?7 and 4 He-7/ bound states considered in this thesis can be searched at 
WASA-at-COSY detection setup in proton-deuteron and deuteron-deuteron fusion re- 
action, respectively. The measurement will be carried out for the beam momentum 
slowly ramped around the r) production threshold corresponding to the range of excess 
energy Q from about -60 MeV to 20 MeV. The existence of the 77-mesic nucleus should 
be visible in the excitation function as a resonance-like structure below the He-?] pro- 
duction threshold. The free //-helium bound states production reactions will be carried 
out in experiment based on measurement of four-momenta of the outgoing particles. 
WASA detector at COSY allows for simultaneous registration of all ejectiles with large 
acceptance, which eg. for the detection of the dd — > ( 4 He-?7)b s — > 3 HepTT~ reaction 
equals about 60%. Spectators from the reactions will be registered mainly in the For- 
ward Detector, while proton-pion pair from the resonance decay will be registered for 
the most part in the Central Detector. The detailed description of the reactions kine- 
matics is given in Chapter 5. 



6.2 COSY-TOF facility 

The measurement of 77-mesic tritium might be realised by means of the quasi-free re- 
action with the time-of-flight spectometer COSY-TOF, a '47r detector' installed at an 
external beamline of the COSY synchrotron. The detection setup is schematically pre- 
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sented in Fig. |6.2| Detailed description of each particular parts and the measurement 
technique might be found in JTTJ [131 1 



Calorimeter 



M ErJangen" 
Startdetector system 




Figure 6.2: Scheme of COSY-TOF detection setup. Charged particles trajectories are reg- 
istered in 'Erlangen' start detector system, while their identification is carried out with the 
Barrel, Ring and Quarrel detectors. The figure is adapted from [13J. 

The proton or deuteron beam accelerated in COSY is extracted and hits a target 
which contains liquid hydrogen or liquid deuterium [Hi [12]. The target is installed 
in front of the 'Erlangen' start detector system consisting of four modular detectors 
(Starttorte, Microstrip detector, Small Hodoscope and Large Hodoscope) designed for a 
precision geometric charged particles tracks reconstruction. The detector covers a polar 
angular range from 3.4° to 74°. The time-of-flight of the charged particles outgoing after 
the interaction of the beam particles with the target, is measured with an accuracy of 
0.25 ns by means of the stop detector situated in cylindric vacuum tank. This detector 
consists of Barrel detector as well as of two Endcap detectors called Quirl and Ring 
and covers angular range from 0.7° to 76.7°. Knowing the time-of-flight and the flight 
length between start and stop detectors, particles velocity is determined. The particle 
momentum can be calculated from the velocity and the mass hypothesis. Additionally, 
neutral particles which take part in the reactions and have not been detected can be 
analysed via momentum and energy conservation. 

The search of ^-tritium bound state can be carried out at COSY-TOF facility 
via the measurement of the excitation function of the nd — > (T-r/)b s — > dpir~ re- 



34 



Experimental setup 



action. The signal from (T-i])b s is expected below the threshold of the nd — > T-n 
production [101 [55]. In the experiment deuteron beam will be used and the nd reaction 
will be analized based on the measurement of the four-momentum of spectator 
proton {psp) [ID] from the dd — > p sp nd —> p sp (T-v)bs — > p sp dpii~ reaction, which was 



schematically shown in Fig 4.2 The advantage of this quasi-free reaction is that the 
Fermi momentum distribution of nucleons inside the deuteron beam allows for the 
scan of energy around the n meson production threshold at a fixed value of the beam 
momentum. Deuterons, protons and pions being products of T-n bound states decays 
will be detected in a multi layer scintillator detectors by measuring their time of flight 
as well as direction. 



7 Simulation results 



In this chapter simulation results of the free and quasi-free r/-mesic bound states pro- 
duction are presented. Monte-Carlo calculations were carried out by means of computer 
programme written in FORTRAN'90 language. 



7.1 Simulation program scheme 

The main purpose of simulations was the determination of the geometrical acceptances 
of WASA-at-COSY and COSY-TOF detectors for free and quasi-free reactions, re- 
spectively and comparing them for different models of nucleon momentum distribution 
inside atomic nuclei and different values of a bound states width. 

In case of free reactions the simulation might be schematicaly described in following 
points for example of reaction (1) which kinematics is presented in details in Chapter 5: 

1. The square of invariant mass of the whole system ^/s^d is distributed randomly 



according to the Breit-Wigner distribution which is given by formula (7.1) and 



shown in Fig. 7.1 




Figure 7.1: Breit-Wigner distribution of square invariant mass -Jsdd 
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N(y/^ = ~ T — 2 (7.1) 



where: m& s = y/sdd thr — B s - a mass of ?7-mesic bound state, 



^Sdd thr = m q + Trujj e - threshold square invariant mass, 



B s - a binding energy of an //-mesic bound state, 
r - width of an 77-mesic bound state. 



In the distribution shown in Fig. 7.1 it is assumed that binding energy equals 
.B s =0.01 GeV and a width is equal to 40 MeV what is in agreement with theoretical 
prediction In simulations a range of V from 7 to 40 MeV was studied. 

2. The N* resonance momentum is distributed isotropically in spherical coordinates 
of ^-mesic nucleus (pp, 9*, <p*) with Fermi momentum distribution of nucleons 
inside 4 He which was presented for three different models in Chapter 3. Next, it is 
transformed into Cartesian coordinates {Pf = {Pf x ,Pp V ,Pp z )) using the following 
equations: 

p** = p *. sinO* -cos(f>* (7.2) 
Pp V = pp ■ sinO* ■ sin<fi* (7.3) 

p* z = p* ■ cos9* (7.4) 

Here the momentum value pp is distributed according to the used model, and the 
direction is simulated isotropically in the space. 

3. The 3 He four momentum vector is calculated (based on spectator model assump- 
tion) in the center of mass frame and transformed using Lorentz transformation 
into laboratory frame. The angle Q^He is also calculated. 

4. Based on JJdd and f>p values, resonance mass m N * is calculated according to 



equation (5.9) 



5. The proton and pion momentum vectors are simulated isotropically in the N* 
frame in spherical coordinates and transformed into Cartesian coordinates. The 



absolute value of p**~ is fixed by equation (5.10). 
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6. The proton and pion four momentum vectors are transformed into the center of 
mass frame and next into laboratory frame by means of Lorentz transformation. 
The angles of outgoing proton and pion in LAB system are calculated. 

7. The histograms of outgoing particles angle distributions and invariant mass dis- 
tribution for all generated events are created. 

8. Knowing angles of outgoing particles and the WASA-at-COSY detector geometry 
it is checked whether generated event can be registered. 

9. The histograms of angular distributions and invariant mass distribution of outgo- 
ing particles for all accepted events are created. 



The simulation of other three free ^-helium bound states production reactions were 
carried out in a similar way. 

The simulation of quasi-free reaction, realised with COSY-TOF detection setup is 
more complex due to the fact that neutron bound in beam deuteron takes part in the 
reaction of ^-mesic tritium formation. The simulation scheme might be presented as 
follows: 



1. The beam neutron momentum is distributed isotropically in spherical coordinates 
with Fermi momentum distribution of nucleons inside beam deuteron (the PARIS 
and CD-Bonn distributions are presented in Chapter 3) and transformed into 
Cartesian coordinates (p*). 

2. The proton spectator four momentum vector as well as its angle with respect to 
the beam direction are calculated in the beam deuteron center of mass frame and 
transformed with Lorentz transformation into laboratory frame. 

3. The neutron four momentum in the beam deuteron frame and in laboratory frame 



is calculated using (4.4), (4.5) and (5.18), (5.19) formulas, respectively. 



The square of bound state invariant mass yjs n d is calculated based on neutron four 
momentum and next events are accepted according to the probability described 
by the Breit-Wigner distribution. 

The next points are analogous like in case of free reaction scheme. Moreover, the 
calculations are carried out for ten values of the deuteron beam momentum f>beam 
in the range from 2.6 GeV/c to 3.5 GeV/c. 



Simulations were conducted for all reactions listed in introduction. Three different 
values of the bound states width: V = {10,25,40} MeV and three different Fermi 
momentum distributions of nucleons inside light atomic nuclei were studied. 
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Simulation results 



7.2 Reactions products - angular distributions 



The ?7-mesic nuclei, which are formed via deuteron-deuteron, proton-deuteron or 
neutron-deuteron fussion, decay according to the mechanism which is described in de- 
tails in Chapter 5. The distribution of momentum vectors of reaction products depend 
on the bound state mass distribution and on the distribution of the nucleons momentum 
inside decayed nuclei. Angular distributions for the outgoing particles being products 
of reaction (1) are presented in Fig. 7.2 
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Figure 7.2: Simulated angular distributions of outgoing 3 He (a), proton (b) and pion (c) 
formed via dd — > ( 4 He-r/)b s —¥ 3 Hep7r~ reaction. Figure shows results for 10 8 generated events 
using the AV18 potential model for the Fermi momentum distribution of nucleons inside 4 He. 



In fact the detectors geometry does not give a possibility to register the whole range 
of outgoing particles angles. An angular ranges covered by respective components of 
WASA-at-COSY detection setup are presented in Fig. |7.2| with shaded areas. 
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7.3 Acceptance 



In order to determine acceptance for respective reaction, the studied excess energy 
range was divided into small intervals. Next for each interval of Q (Q = yfs — y^s thr ) the 
number of events accepted by detector was divided by number of generated events. An 
event is accepted when all outgoing particles being reaction products can be registered 
in the detector. In simulations all possible combinations of particles registration in 
different part of detectors were considered. They are presented for particular reactions 
in Appendix B. 

The acceptance for each reaction was determined for three different values of the 
bound states width T: lOMeV, 25MeV and 40MeV and for three models of nucleons 
Fermi momentum distributions inside 4 He, 3 He or T which are in details described in 
Chapter 3. 

The WASA acceptance as a function of the excess energy Q is presented for two 
reactions of free 4 He-?7 and 3 B.e-r] production: dd — > ( 4 He-r])fe s — > 3 Hep7r~ and pd — > ( 3 He- 
v)bs dpn° — > dp-fy in Fig. 7.3 and Fig. 7.4 respectively. The acceptance is almost 



constant as a function of excess energy and it is independent of the value of T and 
model of Fermi momentum distribution within a few per cent. For two other reactions 
situation is analogous. 
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Figure 7.3: Geometrical acceptances of the WASA-at-COSY detector in case of dd — > ( 4 He- 
rf) hs ->• 3 Hepvr- reaction for the Ti=10 MeV width and AV18 and CDB-2000 models of 
nucleon Fermi momentum distribution inside 4 He (left) as well as for AV18 Fermi momentum 
distribution model and three different values of V (right) . 



40 



Simulation results 



CD 
O 
c 
CO 

-t—l 

Q_ 

8 

o 

CO 



0.6 



0.4 



0.2 















0.6 






anal. form. 


TC6 






r,=10MeV 




CO 


0.4 






AV18 


accept 






r 2 =25MeV 










CDB-2000 






r 3 =40MeV 






0.2 






r 3 =40MeV 

i , . . 




n 





anal. form. 

i , , . 



-60 



-40 -20 

Q [ MeV ] 



20 



-60 -40 -20 

Q [ MeV ] 



20 



Figure 7.4: Geometrical acceptances in case of pd — > ( 3 He-77)& s — > dpir — > dpjj reaction for 
the T3=40 MeV width and three different models of nucleoli Fermi momentum distribution in- 
side 3 He (left) as well as for analitic formula describing nucleon Fermi momentum distribution 
and three different values of T (right). 

The average values of WASA geometrical acceptances for registering the considered 
free reactions were calculated using following formula: 



.4 



N acc {Q) 
Q N gen (Q) 

N 



(7.5) 



where: N acc (Q)- the number of accepted events in a given interval of Q, 



NgeniQ)- the number of generated events in a given interval of Q, 
N- the number of all ranges of summation. 



The obtained results for different gamma values and different models of nucleon mo- 
mentum distributions are presented in following tables: 



[MeV/c] 


AV18 


CDB2000 


n 


0.5296 


0.5380 


T2 


0.5301 


0.5377 


T3 


0.5323 


0.5402 



Table 7.1: Average acceptances of WASA-at-COSY detector for the dd — > ( 4 He-r/)f, s — > 3 Hep7r~ 
reaction for three different T values and AV18 and CDB-2000 models of nucleon momentum 
distribution inside 4 He. 
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[MeV/c] 


AV18 


CDB2000 


n 


0.3785 


0.3889 


T2 


0.3758 


0.3877 


T3 


0.3749 


0.3866 



Table 7.2: Average acceptances of WASA-at-COSY detector for the dd —> ( A ~Re-r])b s — > dppir~ 
reaction for three different T values and AV18 and CDB-2000 models of nucleon momentum 
distribution inside 4 He. 



[MeV/c] 


anal. form. 


AV18 


CDB2000 


n 


0.6175 


0.6071 


0.6139 


T2 


0.6187 


0.6091 


0.6130 


T3 


0.6188 


0.6067 


0.6135 



Table 7.3: Average acceptances of WASA-at-COSY detector for the pd — > ( 3 He-r])bs — > dpir — > 
dp'yy reaction for three different T values and three different models of nucleon momentum 
distribution inside 3 He. 



[MeV/c] 


anal. form. 


AV18 


CDB2000 


n 


0.5054 


0.5059 


0.5052 


T2 


0.5083 


0.5029 


0.5042 


T3 


0.5054 


0.5070 


0.5045 



Table 7.4: Average acceptances of WASA-at-COSY detector for the pd — > ( 3 He-i])b s — > pppir~ 
reaction for three different T values and three different models of nucleon momentum distri- 
bution inside 3 He. 

For the quasi-free reaction dd — > p sp {^-rf)bs PspdpK" the COSY-TOF acceptance 
was calculated for ten of beam momentum values for the range of pb eam =2.6 GeV/c to 
3.5 GeV/c. The excess energy distribution for nd — > (T-r))t, s reactions determined by 



the Fermi momentum distribution of neutron inside deuteron is presented in Fig. |7.5 
for pbeam equal to 2.6 GeV/c, 3.1 GeV/c and 3.5 GeV/c. 

For each case, the number of nd — > (T-n)b s generated events as a function of the 
excess energy Q is represented by solid line while the excess energy distribution for 
accepted events is shown by dashed line. The threshold deuteron beam momentum for 
the dd — > p sp (!T-rj)bs is equal 3.1 GeV/c for the case if the Fermi momentum in the 
beam deuteron is equal to 0. 

Similarly like in case of free reactions, the COSY-TOF acceptance function is ob- 
tained by dividing the number of accepted events by the number of generated events 
for each of excess energy intervals. Respective results for three values of deuteron beam 



momentum are presented in Fig. 7.6 (a), (b) and (c 
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Figure 7.5: The number of generated (solid) and accepted nd — > (T-r])b s — > dpir~ events 
(dashed) as a function of excess energy Q with respect to quasi- free nd — > (T-n)b s reaction at 
deuteron beam momentum pb eam =2. 6 GeV/c (a), pb eam =3.1 GeV/c (b) andp;, eam =3.5 GeV/c 
(c). 



The acceptance of ejectiles registration in the quasi-free reaction is not constant 
function of excess energy. It decreases for Q values corresponding to the Fermi mo- 
mentum of neutron in the deuteron beam equal to zero and results from the fact that 
proton spectator is not accepted by detector geometry if its Fermi momentum value in 
deuteron frame equals p* sp =0. For a comparison, Fig. 7.6 (d) presents that the accep- 



tance dependence of Q would be constant if the proton spectator was accepted. 
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Figure 7.6: Geometrical COSY-TOF acceptances in case of the dd — > Pspi^-fl)bs Pspdpir~ 
reaction at pb eam =2.6 GeV/c (a), pb eam =3.1 GeV/c (b) and Pfe e am=3.5 GeV/c (c). For the 
simulations Paris model of nucleon Fermi momentum distribution inside deuteron and Fi=10 
MeV bound state width and analitic formula describing nucleon Fermi momentum distribution 
inside T were used. Panel (d) represents the reaction acceptance for j>ft ea m=3.1 GeV/c with 
assumption that the detector acceptance for proton spectator registration equals 1. 



In order to determine the effective acceptance for registering the quasi-free 
dd — > p sp (T-r])bs —> p S pdp^~ reaction in the near to threshold Q range from -60 MeV to 
20 MeV when using the TOF detector the following formula is used: 
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Simulation results 



A 



A- 



N, 



(-60,20) 
N gen 



(7.6) 



where: A- the average acceptance given by formula (7.5), 
N gen - the number of generated events, 

-W(-6Q,20) - t ne number of generated events for Qg (-60,20) MeV. 
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Figure 7.7: Effective acceptance for the registration of the quasi- free 
dd — > p sp (T-r])b s — )• p S pdpir~ reaction near the r\ production threshold (Qe(-60,20) MeV) as 
a function of beam momentum. 

The effective acceptance A e rt was calculated for ten values of the beam momentum. 



The dependence A e ff(p beam ) is presented in Fig. 7.7 

In the simulations of the dd — > p sp (T-rj)b s — > p sp dp7i~ reaction it is assumed that the 
square of invariant mass of T-77 nuclei is given by Breit-Wigner distribution. In this 
case an effective acceptance for the registration of the bound state decay products is 
calculated by formula: 



N BW 

ABW _ 4 iV (-60,20) (7 7) 



gen 



where: A- the average acceptance given by formula (7.5) 
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Ng™- the number of generated events accepted with probability 
calculated according to the Breit-Wigner distribution of -Js n d, 



^(-6020)" ^he numrj er of generated events for Qe (-60,20) MeV 
accepted with probability calculated according to the Breit-Wigner 
distribution of J~s^2 



and is presented in Fig. 7.8 as a function of the beam deuteron momentum. 
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Figure 7.8: Effective acceptance for the registration of the quasi free 
dd —> p sp (T-r])b s — )• p S pdpir~ reaction near the rj production threshold (Qe(-60,20) MeV) as 
a function of beam momentum assuming a Breit-Wigner distribution of y/s n d for the bound 
state width r=10MeV. 

From the dependence shown in above figure it results that the highest probabil- 
ity for the registration of considered quasi-free dd — > p S p{pF — 0) (T-77){, s — Y p sp dpiz~ 
reaction is for the beam momentum Pb eam =3.1 GeV/c corresponding to the 
dd — > PspiT-rfbs Pspdp7r~ reaction threshold whereas for the beam momentum above 
and below the threshold the effective acceptance decreases. Above considerations 
were carried out for the Paris model of nucleon Fermi momentum distribution inside 
deuteron, ri=10 MeV bound state width and analitic formula describing nucleon Fermi 
momentum distribution inside T. The calculations are consistent with the one for CD- 
Bonn model of nucleon momentum distribution inside deuteron with an accurancy bet- 
ter than 4% and are independent of models of nucleon Fermi momentum distribution 
inside tritium. 



8 Summary and conclusions 



The main aim of this thesis was to test the feasibility of the measurement of 77-helium 
mesic nucleus production via reactions: 

dd -» ( 4 He-r]) 6s ->■ 3 Hep7r~ 
dd — > ( 4 He-?7)fe s — > dppix~ 
pd — > ( 3 He-r])5 S — » dpn° — > dpjj 
pd ( 3 B.e-i]) bs ->• pppn~, 

with WASA-at-COSY detector and test the feasibility of the 77- T bound state produc- 
tion via reaction: 

nd — > (T-r]) bs — > dpix~ , 

realised through dd — > p sp (T-r])b s — > p sp dpn~ quasi-free reaction with COSY-TOF de- 
tection setup. 

The reaction kinematics for free and quasi-free ?7-mesic bound states production as 
well as the spectator model assumptions were analysed and discussed. Moreover, the 
nucleon momentum distribution inside 4 He, 3 He, T and d, was presented for different 
models based on the analytic formulas and theoretical analysis for each of nuclei. Nu- 
merical data for those distributions are presented in Appendix A. The simulation were 
carried out and acceptances of WASA-at-COSY detection systems for free processes 
were calculated and compared for three assumed values of bound states width as well 
as for three different models of nucleon Fermi momentum distribution. In case of quasi- 
free reaction the effective COSY-TOF acceptance as a function of beam momentum 
was determined. 

The simulation results show that the acceptance as function of excess energy for 
the free reactions of ?7-helium production is in good approximation constant near the 
rj production threshold and is independent of the bound state width value and model 
of Fermi momentum distribution. The calculations present that the absolute value 
of acceptance depends on the reaction channel. The most probable is registering of 
particles outgoing from the pd — > ( 3 He-r])^ — > dpir — > dpy-f reaction. The average value 
of WASA geometrical acceptance is then equal to A=0.61. 

For the quasi- free reaction dd — > p sp (T-Ti)b s — > p sp dpii~ of ^-tritium formation the 
COSY-TOF effective acceptance is dependent on the beam momentum value. It reaches 
the highest value for the beam momentum of j9fc eam =3.1 GeV/c corresponding to the 
dd — > p S p{pF = 0)nd — > p S p(pF = 0)T-^ reaction threshold. Thus the measurement of 
quasi-free reaction products is the most efficient at the 77 production threshold. 

In this thesis the geometrical acceptance for particular reactions was calculated. In 
the future, the efficiency of registration for the outgoing particles in each of reaction 
will be determined. 
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Finally the experiments in which ?7-nuclei will be searched are interesting because 
of following issues: 

1. Potential of discovering of 77-mesic bound states. Observation of that state would 
allow to investigate interactions between the 77 meson and the nucleons inside a 
nuclear matter. 

2. 77-mesic bound systems would provide the information about A r *(1535) [5] reso- 
nance and rj meson properties in nuclear matter [6]. 

3. The existence of the bound states could give a possibility to study a flavour singlet 
component of the quark-gluon wave function of the 77 meson [HI U\. 

The simulations presented in this thesis show that the ?7-helium nuclei measurement 
can be carried out at WASA-at-COSY detection setup, while the 77-tritium bound states 
might be searched at COSY-TOF detector. As a result of simulation it is established 
that the most efficient measurement of dd — > p 3p (T-r))b s — » p sp dp7c~ reaction at COSY- 
TOF detector can be done at the beam momentum of 3.1 GeV/c. 



A Fermi momentum distributions-numerical data 



In Chapter 3 Fermi momentum distributions of nucleons inside d, T, 3 He and 4 He nuclei were 



described and presented in Fig. 3.1, 3.2, 3.3 Numerical data for these distributions are given 
in the following tables: 
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Table A.l: Fermi momentum distribution of nucleon inside deuteron according to Paris model 
and CDBonn model (Chapter 3). The distributions were normalized to unity in the momen- 
tum range from to 400 MeV/c. 
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Table A. 2: Fermi momentum distribution of proton inside 3 He according to three different 
models: analytic formula (equation 3.1) (left), AV18 NN model+UrbanalX three nucleon 
interaction (middle) and CDB2000 NN model+Tucson-Melbourne three nucleon interaction 
(right). The distributions were normalized to unity in the momentum range from to 0.4 
GeV/c. 
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Table A. 3: Fermi momentum distribution of neutron inside T according to three different 
models: analytic formula (equation 3.1) (left), AV18 NN model+UrbanalX three nucleon 
interaction (middle) and CDB2000 NN model+Tucson-Melbourne three nucleon interaction 
(right). The distributions were normalized to unity in the momentum range from to 0.4 
GeV/c. 



B Acceptance-supplement 

The following tables show all possible combinations of particle registration for each of consid- 
ered reactions. Abbreviations 'FD', 'CD-MDC and 'CD-EC denote Forward Detector, Mini 
Drift Chamber and Electromagnetic Calorimeter of WASA-at-COSY detector, respectively. 
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4.3646 


5.6802 


5.6728 


14 


4.6062 


5.4062 


5.4463 


n it: 
U. ±0 


4.7782 


5.1322 


5.1420 


1 £ 
U. lu 


A C TO A 

4.8 /84 


4.7547 


4.8117 


17 
U. ± ( 


4.9075 


4.3501 


4.4151 


18 


4.8685 


3.9455 


4.0184 


1 Q 

u. ±y 


4.7666 


3.5502 


3.6288 


90 

u.zu 


4.6087 


3.1561 


3.2422 


91 


a a non 


2.7o2U 


2.8557 


99 


4.1581 


z.4zoy 


2.5326 


9*i 


3.8835 


2.1320 


0000 


94 


3.5881 


1.8356 


1.9110 


9^ 


0.2oU/ 


1.5385 


1.6683 


9fi 
U.ZO 


2.9092 


1.3591 


1.4551 


97 
U.Z / 


2.00U/ 


1.1797 


1.2418 


98 
U.Zo 


2.3611 


1.0003 


1.0518 


90 

u.zy 


2.0754 


0.8093 


0.9250 


^0 


1.8073 


U. 7 oo7 


n TOO Q 

U. (yoo 


^l 


1.5593 


U.o4oy 


0.6716 


^9 


1.3332 


0.5602 


0.5779 


^"i 
U.oo 


1.1298 


0.4734 


0.5121 




n n a on 


0.4194 


0.4463 




U.OO 


0.7902 


r\ or 1 

U.3868 


0.3805 


^fi 
U.OU 


0.6523 


0.3542 


VJ.OOOZ 


^7 


0.5339 


0.3216 


U.oU2o 


^8 
U.OO 


0.4333 


U.2889 


U.2 (2[) 


^0 


0.3488 


0.2759 


0.2415 


0.40 


0.2784 


0.2665 


0.2149 


0.41 


0.2204 


0.2570 


0.2012 


0.42 


0.1730 


0.2475 


0.1875 


0.43 


0.1348 


0.2380 


0.1737 


0.44 


0.1041 


0.2331 


0.1601 


0.45 


0.0798 


0.2295 


0.1510 


0.46 


0.0607 


0.2260 


0.1439 


0.47 


0.0458 


0.2224 


0.1368 


0.48 


0.0343 


0.2189 


0.1297 


0.49 


0.0254 


0.2153 


0.1226 


0.50 


0.0187 


0.2133 


0.1177 



Table A. 4: Fermi momentum distribution of neutron inside 4 He according to three different 
models: analytic formula (equation 3.2) (left), AV18 NN model+UrbanalX three nucleon 
interaction (middle) and CDB2000 NN model+Tucson-Melbourne three nucleon interaction 
(right). The distributions were normalized to unity in the momentum range from to 0.5 
GeV/c. 
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particle 


case 1 


case 2 


case 3 


case 4 


3 He 


FD 


FD 


FD 


FD 


proton 


FD 


CD-MDC 


FD 


CD-MDC 


7T~ 


FD 


CD-MDC 


CD-MDC 


FD 



Table B.l: Possibilities of particles registration in the dd — > ( 4 He-7?){, s — > 3 Hep7r reaction. 



particle 


case 1 


case 2 


case 3 


case 4 


case 5 


deuteron 


FD 


CD-MDC 


FD 


FD 


FD 


protonl 


FD 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


proton2 


FD 


CD-MDC 


CD-MDC 


FD 


CD-MDC 


7T~ 


FD 


CD-MDC 


CD-MDC 


FD 


FD 



case 6 


case 7 


case 8 


case 9 


case 10 


case 11 


FD 


FD 


FD 


FD 


CD-MDC 


CD-MDC 


CD-MDC 


FD 


FD 


FD 


CD-MDC 


FD 


FD 


CD-MDC 


FD 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


FD 


FD 


CD-MDC 



case 12 


case 13 


case 14 


case 15 


case 16 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


FD 


FD 


FD 


CD-MDC 


FD 


FD 


CD-MDC 


FD 


FD 


CD-MDC 


CD-MDC 


FD 


FD 


FD 



Table B.2: Possibilities of particles registration in the dd — > ( 4 He-r/)f, s — > dppir reaction. 



particle 


case 1 


case 2 


deuteron 


FD 


FD 


proton 


CD-MDC 


FD 


gammal 


CD-EC 


CD-EC 


gamma2 


CD-EC 


CD-EC 



Table B.3: Possibilities of particles registration in the pd — > ( 3 He-r?)f, s — > djyyy reaction. 
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particle 


case 1 


case 2 


case 3 


case 4 


case 5 


protonl 


FD 


CD-MDC 


FD 


FD 


FD 


proton2 


FD 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


proton3 


FD 


CD-MDC 


CD-MDC 


FD 


CD-MDC 




FD 


CD-MDC 


CD-MDC 


FD 


FD 



case 6 


case 7 


case 8 


case 9 


case 10 


case 11 


FD 


FD 


FD 


FD 


CD-MDC 


CD-MDC 


CD-MDC 


FD 


FD 


FD 


CD-MDC 


FD 


FD 


CD-MDC 


FD 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


FD 


FD 


CD-MDC 



case 12 


case 13 


case 14 


case 15 


case 16 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


CD-MDC 


FD 


FD 


FD 


CD-MDC 


FD 


FD 


CD-MDC 


FD 


FD 


CD-MDC 


CD-MDC 


FD 


FD 


FD 



Table B.4: Possibilities of particles registration in the pd — > ( 3 He-r/)f, s — > pppir reaction. 



particle 


case 1 


proton S p 
deuteron 
proton 

7T~ 


START and STOP DETECTORS 
START and STOP DETECTORS 
START and STOP DETECTORS 
START and STOP DETECTORS 



Table B.5: Possibilities of particles registration in the dd — > p sp (T-r/)fe s — > p sp dpir reaction. 



